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We study the stability of the spin gap phase in the (7(1) slave-boson theory of the t — J model in 
connection to the underdoped cuprates. We approach the spin gap phase from the superconducting 
state and consider the quantum phase transition of the slave-bosons at zero temperature by intro- 
ducing vortices in the boson superfluid. At finite temperatures, the properties of the bosons are 
different from those in the strange metal phase and lead to modified gauge field fluctuations. As a 
result, the spin gap phase can be stabilized in the quantum critical (QC) and quantum disordered 
(QD) regime of the boson system. We also show that the regime of QD bosons with the paired 
fcrmions can be regarded as a strong coupling version of the recently proposed nodal liquid theory. 
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The pseudogap behavior of underdoped cuprates in 
various physical properties has been a subject of inten- 
sive research recently H. The importance of the subject 
comes from the fact that understanding this behavior 
and relating it to the properties of the superconducting 
state may give us an important clue for the mechanism 
of the superconductivity. Among many different theories 
for the pseudogap, one of the earliest proposals was the 
mean field theory of the t — J model which is an effective 
low energy theory of the Hubbard model in the limit of 
large on-site Coulomb repulsion 0|. 

In order to explain why this mean field theory is ap- 
pealing, let us begin with the slave boson representation 
of the t — J model. In the strong coupling limit, the 
double occupancy of the electrons at each site is prohib- 
ited, thus it is convenient to describe the Hilbert space 
of the electrons in terms of a neutral spin 1/2 fermion, 
Ha ' representing the singly occupied sites with a spin up 
or spin down electron, and a spinless charge e boson, bi, 
keeping track of empty sites. As a result, the electron op- 
erator can be written as cj = fj bi with the constraint 



Eo fLhc, + b\bi = 1. The t — J model can be written as 
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where = f 3 A ■ dl and Ao represent the external vec- 
tor and scalar potentials. Here aoi is the Lagrange multi- 
plier enforcing the local constraint. Using the Hubbard- 
Stratonovich transformation, the action can be rewritten 
as 
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In the mean field theory, it was 



found that \Qij \ = Q = const, and a,j =0. It has been 
established that there are four different phases at the 
mean field level depending on the values of Ay and (bi) 
§. i) Superconducting phase: A„ = (e^fjjjp) + 



and (bi) 7^ 0. ii) Spin gap phase: Aij ^ and (bi 



0. 



iii) 'Fermi liquid' phase: Aij = and (bi) ^ 0. iv) 
'Strange metal' phase: Ay = and (bi) = 0. As the elec- 
tron is a combination of the fermion and boson, an excita- 
tion gap of the electron will be generated from the gap of 
the spin carrying fermions in the spin gap phase. There- 
fore, this theory already suggests the spin gap phase as a 
possible candidate for the pseudogap behavior of under- 
doped cuprates. 

However, it was later found by Ubbens and Lee that 
the spin gap phase of the mean field theory is unstable 
against the fluctuations about the mean field state when 
the spin gap phase is approached from the strange metal 
phase [|| . One can see from Eqj^ that the fluctuation in 
the phase of Qij can be represented in terms of a U(l) 
gauge field a through ay = J? a • dl. This is associ- 
ated with the internal U(l) symmetry in the slave-boson 
representation of the t—J model. The previous study 
mentioned above would imply that the U(l) slave-boson 
theory does not support the spin gap phase. 

In order to make the later discussion more concrete, let 
us first reproduce the arguments of Ubbens and Lee [[| . 
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Upon approaching the spin gap phase from the strange 
metal phase, one can evaluate the free energy cost for 
opening up the gap for the fermions in the following fash- 
ion. First we evaluate the contribution to the free energy 
from the gauge field fluctuations in the strange metal 
phase. 
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where -D(q, lo) is the gauge field propagator and given 
by D-^q.w) = Ilf(q,a;) +n^'(q,w). Here H£? and 
II^ are the fermion and boson current-current correla- 
tion functions respectively. In the strange metal phase, 
Hp and H J g were assumed to have the free fermion 
and boson forms: U JJ — 
XbQ 2 - Here 7 = 2n e /k 

^Tbe/t i)/(247rms), where n e is the electron density, 
nip = \/(2JQ), tub = l/(2iQ), and Tbe = 2irx/mB 
(x is the doping concentration) is the mean field boson 
condensation temperature. As a result, D~ l is given by 
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where x — Xf + Xb- Using EqJ| and Eq.|[ we get 
F g oc T 5 / 3 . When the gap is opened, A can be intro- 
duced to cutoff the frequency integral and the effect is 
simply replacing T by A in F„. Thus, the free energy 
cost for opening up the spin gap is proportional to A 5 / 3 . 
Since the mean field pairing energy gain is proportional 
to — A 2 , the free energy cost from the gauge field always 
dominates. Thus the spin gap phase cannot be stabilized. 
In their work, the bosons in the strange metal phase were 
assumed to behave rather classically. The transition from 
the strange metal phase to the spin gap phase occurs due 
to the pairing of the fermions while the bosons were as- 
sumed to be still classical. 

In this paper, we suggest that the spin gap phase can 
be stabilized at low temperatures because the properties 
of the bosons are different from those in the high tem- 
perature strange metal phase. Motivated by the phase 
diagram of the cuprates, we suggest that at zero temper- 
ature there exists a quantum disordering phase transi- 
tion of the bosons driven by vortices while the fermions 
remain paired across the transition (Fig. 1). The low 
density boson system (x < x c ) becomes an insulator due 
to the condensation of vortices. On the other hand, if 
x > x Cl the bosons are in the superfluid state while the 
vortices are in the insulating state. Recall that the elec- 
tron resistivity is given by p e — pp + ps, where pp and 
Pb are the fermion and boson resistivities respectively. 
Since pp = and ps diverges when x < x c , the zero tem- 
perature spin gap phase is an insulator. We found that, 
near the QC point of the phase transition, the proper- 
ties of the bosons are significantly modified. Using the 



boson-vortex duality, we obtain the current-current cor- 
relation function, H^, for the QD bosons. Due to the 
modified bosonic properties, the gauge field propagator 
has a different form at low temperatures, which leads 
to, according to Eq.||, the free energy cost for opening 
the spin gap smaller than the pairing energy gain, —A 2 . 
Thus the spin gap phase can be stabilized. In addition, 
we suggest that the recently proposed nodal liquid theory 
can be regarded as a weak coupling version of our theory 
if the constraints are ignored B . 

Let us begin with the superconducting state where the 
fermions are paired (Ay ^ 0) and the bosons are in the 
superfluid state ((hi) ^ 0). The phase diagram of the 
cuprates tells us that the superconducting state exists 
when the doping concentration is larger than a particu- 
lar value x c . The question is whether the normal state 
in the regime x < x c can be obtained by suppressing 
the boson superfluid while the fermions are still paired. 
In order to answer this question, we have to understand 
the nature of the quantum disordering phase transition 
of the bosons across x c . We suggest that the latter is due 
to the condensation of vortices. The vortices carry the 
flux quantum hc/e because the bosons carry charge e. In 
order to describe the QD state of bosons, it is convenient 
to use the dual representation of the bosons || . Notice 
that the external electromagnetic fields only couple to 
the bosons. The continuum limit of Cb in the presence 
of the external fields is given by 

C B = bUd T - ia - iA )b - -^—bUw iA) 2 b . 
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Following Ref. ||, let b = yfpcj) where p is positive defi- 
nite, which corresponds to the boson density, and <j> is a 
unimodular complex field satisfying <j>*(j) = 1. Then the 
action (Eq.||) becomes 
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Next, we decouple the second term by a Hubbard- 
Stratonovich transformation, which leads to 
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where J M = (p, J) is the boson three-current and = 
(A ,A). In order to isolate the vortices, we write (f> as 
4> = (f> v e l6 , where 6 is single valued and (f> v represents 
the vortices. Integration over 9 gives the continuity con- 
straint d^J^ — 0. One can solve this constraint by in- 
troducing a new field M\ through J M = e^ v \d v M\ = 
(d x M) M . It is also useful to introduce SM^ = M M — M M 
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with (d x M)o = x and (d x M)i,2 = 0, where x is 
the average boson density. Then the action in the long 
wavelength limit can be written as 



C B = ^(\(d x SM)^ 2 + |(9 x 5M) 2 \ 2 ) 
+iSM^{J» -(dx a)„ - (8 x A)fj] 
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where (9 x a)^ = e^uxd^ax and is the vortex three- 
current = e^xdv^l^^v = (d x 4>lj(l> v )v When 
the vortices are condensed, the vortex transverse current- 
current correlation function in the long wavelength and 
low energy limit should have the form: (J v J") — 
Cp s v (T, x), where p s v {T, x) is the superfluid density of the 
vortex condensate and C is a constant. Using this cor- 
relation function and integrating out M M degrees of free- 
dom, the effective action for the transverse part of a^+A^ 
field becomes 
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Comparing with Eq||, we can read off the boson current- 
current correlation function at finite temperatures: 
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Since p s v becomes smaller and smaller as the critical point 
is approached, if we write down II^ as — XbQ 2 with 
Xb fx we have \f ^ Xb near the critical point. 

Then the gauge field propagator at finite temperatures 
can be written as 
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Using the above propagator and Eq.[|, we obtain 
F s cx[p*(7>)] 2 / 3 T 5 / 3 . 



(11) 



(12) 



Now some remarks on the temperature dependence of 
p%(T,x) are in order. Using dimensional analysis of hy- 
perscaling, i.e. the free energy within a correlation vol- 
ume £ d £ z is non-singular in the vicinity of the quantum 
critical point, the scaling form of p s v {T, x) in d-dimensions 
can be written as H 



pl{T,x)=e- d - z F((, z T) , 
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where £ oc \x — x c \~ v is the correlation length of the vor- 
tex superfluid and F(x) is a scaling function. Here z 
is the dynamical critical exponent and v is the correla- 
tion length exponent. If £ 2 T 3> 1, i.e. in the quantum 
critical regime, the vortex superfluid density must be in- 
dependent of £ and only depend on temperature. In two 
dimensions, this implies F(x) ~ x for x ^> 1, leading to 
pi oc T. Substituting the latter into Eq.|l2|, we find that 
the free energy cost for opening up the spin gap due to 



the gauge field becomes F g ~ T 1 l z in the QC regime. 
Similarly, in the quantum disordered regime, £ Z T <C 1, 
the scaling function F(x) — > const for x <§; 1, p% becomes 
temperature independent p% oc '. In this case, Eq.[T^ 
becomes F g oc £~ 2z /3y5/3^ wh en T < A, one would use 
A as a low frequency cutoff instead of T. As long as 
A > B\x — x c \ zv , one can trade T and £ -z with A in the 
above expression to obtain the energy cost for opening up 
the gap F g oc A 7 ' 3 . Thus the mean field fermion pair- 
ing energy gain, —A 2 , wins and the spin gap phase can 
be stabilized in both the QC and the QD regimes. We 
stress that this is a qualitatively very different behavior 
compared to the previous claim pi. 
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FIG. 1. Schematics of the phase diagram. The dashed 
lines indicate the crossover between QC and QD, QC and 
RC, and QC/RC and the strange metal (SM) regimes. 
The spin gap is stable in the QC and QD regimes. 

At sufficiently high temperatures, the bosons would be- 
have classically. In this regime, the argument of Ubbens 
and Lee can be applied Therefore, there should be 
a crossover from either QD or QC regime to the high 
temperature strange metal phase where the gauge field 
fluctuations destroy the spin gap phase. Fig. 1 shows a 
schematic phase diagram in the T — x plane. Notice that 
there exists a crossover line between the strange metal 
phase and the spin gap phase stabilized in both the QD or 
QC regime of bosons. The dashed lines (T = B\x — x,\ zv 
line) represent the crossovers between the QD, QC, and 
renormalized classical (RC) regimes of bosons. Since the 
bosons in the RC regime behave essentially in a classical 
fashion, the spin gap phase is not stable there. Thus, the 
spin gap phase is restricted to the underdoped regime. 
The crossover from the QC regime to RC regime occurs 
when T < B(x — x c ) zu and x > x c . Thus the transi- 
tion from the strange metal phase to the superconducting 
phase is likely to occur in the RC regime (i.e., T c may 
be lower than B(x — x c ) zv ) so that it is characterized by 
a more or less classical transition of the bosons to the 
superfluid state. On the other hand, if T < B(x c — x) zv 
and x < x c , the crossover from the QC regime to the QD 
state of bosons occurs. Thus the antiferromagnetic phase 
is in the quantum disordered regime of the bosons. In the 



3 



low temperature limit, the relevant excitations in d-wave 
spin gap phase would be the neutral Dirac fermions. At 
lower doping concentrations, the SU(2) fluctuations |t]] 
may become important and these neutral Dirac fermions 
may be confined when x is sufficiently small and the an- 
tiferromagnetism occurs JjJU ■ 

Now let us discuss the relation between the present the- 
ory and the nodal liquid theory proposed recently, where 
the quantum disordering of Cooper pairs was discussed 
[pj . Let us begin with the continuum limit of Cf assum- 
ing that Aij — A is a uniform complex number and we 
will comment on the d-wave case later: 
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Let A = |A|e^. After doing the gauge transformation 
fa — * f a e 1 ^^ 2 and b — ► 6e 4 ^/ 2 , we obtain the action in the 
absence of the external fields: 



C 



fl(d T - (i - z~a )f a - 2^4(V - i*?f a 
+|A|(e^/ a + e^/t/tf, 

C B = bUd T - ia )b - —— & f (V - ia) 2 b , 
zrriB 



(15) 



where 5,^ = a M + (<3„0)/2. This amounts to fixing the 
gauge for et M . In the superconducting state, the order 
parameter in the strong coupling limit can be written as 



( £Q/3 I n the gauge choice we 



have taken, (e a/3 f a fp) is always real. Therefore, in the 
superconducting state, we get A e w |A|(66) « |A|6ge l9 . 
Here 6 = boe lB / 2 is taken. Notice that the phase of A is 
dictated by the phase of the bosons and |A e | = \A\b$. In 
the superconducting state, we substitute b by b^e 1 ®! 2 in 
the action. Integrating out ao,a, we obtain 
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where = (p F ,J F ) with ^ = and J F = 

2^(f a Vf a - H.c). n£„ = (J,f J^) is the fermion 
three-current correlation function. Here g — 1/2 and 
gi 2 = ( + Ilz*l]. Also k = A^(O)- 1 = an d 

Ki : 2 = ^(^r^").9i,2- The continuum action derived above 
can be also obtained directly from the lattice model given 
by Eq.0§. 

If A^ has the d-wave symmetry, the relevant degrees of 
freedom in the low energy limit are the excitations near 



the nodes. These excitations have the Dirac spectra and 
can be represented by d a = (d a -[, d a \), where a = 1, 2 rep- 
resent two pairs of nodes. Following Ref. jl, the fermion 
part of the action becomes 

Cd = d\(d T — i'FT z id x — v&T x id y )d\ 

+d\{d T - v A T z id y - v F T x id x )d 2 , (17) 

where t x ' v ' z are the Pauli matrices and va = A/y/2. The 
coupling between the fermion and the phase fields is also 
changed to Cd,e — 9^Jt{dnO), where j£ is the current 
of the Dirac fermion d. In this case, the corresponding 
action in Eq.|l^ looks similar to that of the nodal liq- 
uid theory if the additional terms coming from the con- 
straints are ignored. Thus we conclude that our theory 
can be also interpreted as the strong coupling version of 
the nodal liquid theory. However, due to the presence 
of the constraints, the structure of the theory is not ex- 
actly the same. Notice also that the presence of the hc/e 
vortices as well as the hc/2e vortices was pointed out by 
Sachdev |l0| as a consequence of the constraint imposed 
by the gauge field. We believe that the presence of the 
hc/e vortices suggested by the nodal liquid theory comes 
out naturally in the present strong coupling theory. It 
was also suggested that the vortex condensate supports 
charge e solitons (dual vortex of the vortex condensate) , 
holons in the dual picture. Using the boson-vortex 
duality, one can see that the slave-bosons in the boson 
picture may be the natural candidates for the solitonic 
excitations in the dual representation. 

In summary, using the t/(l) slave-boson representa- 
tion of the t — J model, we obtained the spin gap phase 
by quantum disordering the slave-boson superfluid of the 
superconducting phase while the spin-carrying neutral 
fermions are paired. We found that the spin gap phase 
can be stabilized in the QC and QD regime of the bosons 
against the fluctuations about the mean field state. We 
also showed that the spin gap phase at zero tempera- 
ture obtained in this fashion can be regarded as a strong 
coupling version of the nodal liquid phase. 
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